Optical transparency of graphene as determined by the fine-structure constant 
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The observed 97.7% optical transparency of graphene has been linked to the value 1/137 of the 
fine structure constant, by using results for noninteracting Dirac fermions. The agreement in three 
significant figures requires an explanation for the apparent unimportance of the Coulomb interaction. 
Using arguments based on Ward identities, the leading corrections to the optical conductivity due 
to the Coulomb interactions are correctly computed (resolving a theoretical dispute) and shown to 
amount to only 1-2%, corresponding to 0.03-0.04% in the transparency. 



The optical transparency of graphene is determined by 
its optical conductivity a {uj) and c, the speed of lighiji: 



t{uj) = (l + 27r(T (cj)/c)" 



(1) 



Recent experimentsS^ on suspended graphene found 
t {uj) ~ 0.977, independent of w, in the visual regime 
(450nm < A < 750nm). This observation (see also 
Refs. can be elegantly rationalized in terms of 

non-interacting Dirac particles with optical conductiv- 
ity£ cr(°) (fcsT < w < D) = f eV/i • Here D is the upper 
cut off energy for the linear dispersion, of order of sev- 
eral electron volts and T is the temperature. Assuming 
cr = crC^) yields t{u)) ~ 0.9774629(2), in excellent agree- 
ment with experiment. Thus, the optical transparency 
of non-interacting graphene t (w) = (1 -I- ttqiqed/S) is 
solely determined by the value of the fine structure con- 
stant of quantum electrodynamics: q;qed = e^/ {he) ~ 
1/137.035999(6). Despite the beauty of this reason- 
ing, a natural question emerges: Why can one ignore 
the electron- electron Coulomb interaction? After all 
the Coulomb interaction in graphene is poorly screened 
and its strength is governed by its own, effective fine 
structure constant a = e^/ (hv) ~ 2.2 that is signifi- 
cantly larger than q;qed because of the smaller velocity^ 
V ~ lO^m/s. The quantitative agreement between ex- 
periment and a non-interacting theory clearly requires a 
quantitative analysis of the size of interaction corrections 
to the optical conductivity and transparency of graphene. 

In this Brief Report we determine the leading interac- 
tion corrections to the optical transparency and demon- 
strate that they amount to only 0.037% in the visual 
regime. This surprisingly small correction is the conse- 
quence of i) a perfect cancellation of the divergent (i.e. 
proportional to lnD/|ti;|) parts of Feynman diagrams 
that contribute to the conductivity and ii) a near cancel- 
lation of the non-divergent contributions. While the first 
result has been stated earlier by us^ as well as in Ref . i, 
the latter effect has been a subject of a dispute^ i^°i^^ . 
Below we resolve this dispute and demonstrate that the 
leading perturbative correction to the conductivity was 
correctly analyzed by Mishchenko in Ref. [l^. We show 
that perturbative corrections to the conductivity must be 
obtained by guaranteeing that momentum cutoffs, used 
to regularize divergences, are introduced in a fashion that 
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Figure 1: (Color Online) Optical transparency, Eq. ([T]), of 
graphene, from Ref. Q (points) along with theoretical curves 
for the case of interacting graphene within the present the- 
ory, Eq. ^ with Ci ~ 0.01 (solid red line), according to the 
theory of Ref. iS (Ci ~ 0.51; dot-dashed green line), and for 
noninteracting Dirac fermions (dashed blue line). 

respects Ward identities and thus guarantees charge con- 
servation. 

The low energy Hamiltonian for electrons in 
graphenelS is obtained by expanding to leading order in 
gradients near the nodes of the tight-binding dispersion, 
yielding the following nodal-fermion Hamiltonian: 



p{r)p{r') 



(2) 



Here i = 1 • • • A'^, with TV = 4 counting the two spin in- 
dices and two independent nodes in the Brillouin zone, cr^ 
are the Pauli matrices, and we have set /i = 1. It is simple 
to show that the charge density p (r) = X]t=i ^1 (i") V'i (r) 
and current density j (r) =i'X]i=i V'J (r) (ypi (r) are re- 
lated by the continuity equation: 



dp 

dt 



V-j = 0. 



(3) 



The optical conductivity is related by the Kubo formula 
to the retarded current-current correlation function: 



cr(w) 



2uj 



Im[/f,(q = 0,c^)-t-4^^(q = 0,L^) 



(4) 
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Here, fj^^,{c\,bj) is the retarded correlation function 
determined from the Matsubara function f^v{Q) = 
(j/i {Q) jv {—Q)) by analytically continuing iO. u+iO^ . 
We use the convention Q — {—i^, q) and correspondingly 
write jo {Q) — P [Q) for the charge density. 

The theory for the optical conductivity in graphene 
with electron-electron Coulomb interaction was devel- 
oped in Refs. [SiS- Using renormalization group (RG) ar- 
guments it holds that the effective fine-structure constant 
of graphene, a, becomes a running coupling constant a (l) 
where I is the flow variable of the RG approach. In the 
case of graphene a{l) decreases logarithmically as one 
lowers the typical energy scale^'^'i^'i^ii^'i^iii. The opti- 
cal conductivity a{uj,T,a) at frequency w, temperature 
T and for the physical coupling constant a is related to its 
value at a rescaled frequency lor (l) = Z (/) u, rescaled 
temperature (l) — Z (l) T as well as the running 
coupling constant via 



a (w, T, a) = (7 {ujR (l) , Tr (/) , a (l)) . 



(5) 



The scaling factor up to one loop is given by Z (l) — 
e~' (l -I- jl). Equation ^ implies that the conductiv- 
ity is scale invariant under the RG flow. This result is 
true to arbitrary order in perturbation theory as can be 
shown following arguments by Grossi^. It is physically 
due to the fact that the electron charge is conserved^-. 
The scaHng functions ujr (l) and Tr (l) grow under renor- 
malization while a (l) decrease a^i^^ . Thus, in the relevant 
collisionless regime w ^ T is it sufficient to analyze the 
high frequency {lur (l) ~ D), weak coupHng limit where 



a{D,0,a) 



[l+Cia + C2a^ + ■■■] 



(6) 



Here, the numerical coefficients Ci are determined by per- 
forming an explicit perturbation theory. The scaling law 
Eq. yields the conductivity as a function of frequency 
where we replace a by the running coupling constant 



(7) 



here 



obtained 
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to leading logarithmic 
accurac y"'"'""'""'""'""'"' . The result is that interac- 
tions only give rise to additive corrections to cr^"^ that 
are of the form 



cr (w) = cr(") [1 + Cia (uj) + Cja^ (lu) + 



(8) 



Note, this behavior is correct in the collisionless regime 
^ ^ ksT. Qualitatively different behavior occurs in the 
opposite, hydrodynamic regimei^ uj <^ ksT. 

Since a (lj — > 0) = 0, it follows from Eq. ^ that 
tj (a; 0) (t(°) . However, a {uj) only vanishes as 
4/log(D/w) and corrections could easily be significant 
in the visible part of the spectrum where u and D are 
comparable. The dominant correction is due to the Cia 
term and will be analyzed in this Brief Report. Com- 
bining Eqs. ^ and |[8]), and neglecting the higher-order 
terms, we have 



cr(w) 



1 



Cia 



l + \a\og{D/u) 



(9) 
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Figure 2: Feynman diagrams for the leading order contribu- 
tions to fij,v{Q) and <j{uj). Diagram A is the 0(0:") contribu- 
tion, while diagrams B, C, D are 0{a). Full lines represent 
fermions and dashed lines represent the Coulomb interaction. 

Calculations of Ci were presented in Refs. Isl and [lol. how- 
ever with different results. While the authors of Ref. 
obtained Ci = (25 - Gtt) /12 ~ 0.513, MishchenkciQ ob- 
tained a significantly smaller value Ci = (19 — 67r) /12 ~ 
0.0125 , which was however disputed in Ref. |Tl|. Deter- 
mining the correct value of Ci is important for two rea- 
sons. First, there is no obvious mistake in either Ref. |3 
or RefllO. It is clearly important from a purely theoret- 
ical point of view to settle this issue and set the criteria 
for correct calculations of interaction effects in graphene. 
Second, as we discuss in more detail below and illustrate 
in Fig. [ll the coefficient determined in Ref. d is not con- 
sistent with experiment, impling that qualitatively new 
phenomena or even higher order corrections would have 
to be invoked to understand the observations of Ref. ^2. 

We now obtain (j[lo) by calculating the correlation 
function ff_iu[Q) which, as follows from Eq. satisfies 



.U. iQ) = 0, 



(10) 



with the repeated index summed over fi — 0,x,y. When 
calculating f^v{Q)^ the leading contributions to which 
are shown in Fig. [2l we must ensure that Eq. (fTO|) is 
satisfied at each order in a. 

The zeroth-order contribution, Fig.[2l\, corresponds to 
the current-current correlation function f^j^) [Q) of non- 
interacting Dirac particles and yieldsi^ 



4°; (Q) 



N 



16IQI 



(11) 



which obeys Eq. (|T0|) . Performing the analytic continua- 
tion and inserting the result into the Kubo formula yields, 
after restoring proper units, a^^^ = ^ for iV = 4. This 
leads to the 97.7% optical transmission discussed above. 

Next we analyze the three leading corrections to ct'"^ 
as shown in Fig. [2l3-D. The diagrams in Fig. and C 
yield the same contribution with interactions entering via 



3 



self energy insertions with the leading self energy 

S(p) = -e'|^/(p-p')G(P'), 

where ... = T^.. J,, ffp- and ^(p) 
Fourier-transformed Coulomb interaction, 
propagator is given by 



(12) 



= § is the 
The fermion 



G{P) 



(13) 



The self energy, Eq. (|T2l) . diverges logarithmically and 
must be regularized, for example by introducing an upper 
momentum cut off A ~ D/v. We will show that the 
discrepancy between previous calculations of (7(0;) can be 
traced to the fact that, in Eq. lfT2l) . there are two obvious 
ways to introduce the ultraviolet (UV) cut-off A. 

Thus, upon evaluating the frequency summation and 
momentum integrals, we obtain 



E(p) = -azjp • cr In 



4Ac 
P 



(14) 



where the number c depends on the cutoff procedure, 
with c = e~^/^ if we evaluate the momentum integral us- 
ing the cutoff Ip'I < A, i.e. by confining fermion states 
to a circle near the node. (Note we always discard con- 
tributions that vanish for K/p — > 00). On the other 
hand, if we evaluate the momentum integral by restrict- 
ing |p — p'l < A, i.e. by having a finite Coulomb inter- 
action at short distances (see also Ref. [ll|) , we find that 
c = e^/^. This corresponds to replacing the Coulomb 
potential via y(p) ^ 1/a(p) = 2tt0{A - |p|)/|p|. We 
emphasize that the log-divergent contribution to S(p) 
is independent of the regularization procedure, with the 
difference being in the subleading contributions. 

As we will show, the two different values that have been 
determined for Ci in Refs. |3 and !10:, respectively, are di- 
rectly related to the two different values for c in the self 
energy as it enters in the diagrams of Fig. [2B and C. 
The diagram Fig. [JjD is unaffected by the regularization 
procedure. Which result for c, i.e. which regularization 
procedure, is correct? The answer comes from the Ward 
identity, which, as we show next, is only satisfied if we 
implement the momentum cutoff by restricting the mo- 
menta in the Coulomb potential to |p — p'| < A, implying 
that Mishchenko's resultiS is correct. 

To demonstrate that the proper cut off procedure is 
to restrict |p — p'| < A in the Coulomb potential, we 
analyze the leading interaction corrections (Fig. [2lB,C, 
and D), which we call f^n) {Q). These satisfy: 

Q^fi'JiQ) ^ NaJ J Va{p-p') (15) 

xTr {G(P' + Q)G{P + Q)G{P' + Q)a, 

-G{P')a,G{P')G{P)} ■ 

This result was obtained from the three diagrams 
Fig. by simply using the identity 

GiP){inao-(i-a)GiP + Q) = GiP)~GiP + Q), (16) 



and the cyclic property of the trace. At this point, the 
UV cutoff only enters via Va(p), so that we can shift 
P ^ P — Q and P' ^ P' — Q in the first term, again use 
the cyclic property of the trace, and obtain 



fil\Q) = o, 



ifj.j p.. 



(17) 



the required result. Note that other regulation schemes 
for the UV behavior will not necessarily work in this 
way. In particular, regulating the momenta by restrict- 
ing the Green-function momentum arguments amounts 
to replacing G{P) G{P)Q{A — |p|); with such a re- 
placement, Eq. ifTG)) and thus Eq. (fTSl will not be valid. 
We conclude, then, that in graphene momenta must be 
regularized using Va(p). 

The same conclusion can be arrived at by considering 
the leading corrections to the current vertex 

A^(P, Q)^~a f G{K)a^G{K + Q)Vk{\<l - p). (18) 

JK 

Again using Eq. (fTBl) . we obtain 



QpA,. = -a I (G{K + Q)-G(K))VA{'k-p) 

JK 



S(P + Q)-S(P), 



(19) 



which is the correct Ward identity. Once again, alter- 
nate schemes for cutting off the momentum integrals are 
not guaranteed to yield a proper Ward identity of this 
form. Our finding that a regularization in terms of a hard 
fermion cut-off violates charge conservation is analogous 
to the observation in QED that incorrect regularization 
schemes yield unphysical results such as a photon mass^S. 

Our final tasks are to evaluate the contributions to 
the current-current correlation function and to determine 
the conductivity using the Kubo formula Eq. (|4]). We 
first consider the diagrams B and C, which are identical. 
Recognizing the self-energy insertion, we have (with an 
overall 2 for the two diagrams): 

/« \bc = -2N jTT[GiP)a,G{P + 0)a.G(P)S(p)]. (20) 

Evaluating the trace, and performing the frequency inte- 
gral and analytical continuation yields 

T .CUB/ M Nauj , 8uAce"^/^ , , 

Im/«« q = 0,c^ BC = -^In , 21 

'^'^ 16 uj 

where the repeated ^ index refers to the sum over the xx 
and yy components as in Eq. (|4]) . Analyzing the diagram 
D of Fig. [21 which can be written as 

fj^'JlD - ^NjTT[G{P)A^{P,Q)GiP + Q)a,l (22) 

it turns out that the result does not depend on the details 
of the regularization procedure and yields 

,1-1 -iR/ M Nauj /, 8vA 19 — 67r\ , , 

Im/W^(q = 0,^)b ^ ( 1'^ — + ■ (2^) 
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By examining Eqs. (|2T1) and (|23l) . it is clear that the de- 
pendence on the high energy scale A vanishes, in agree- 
ment with general scaHng arguments^>^. Plugging these 
results into Eq. |[4]) yields Eq. ^ with coefficient 



C,^l^-ilnce~V^ 
12 2 



(24) 



We indeed see that the correct cut-ofF procedure, with 
c = e^/^, yields Ci = , whereas the other cutoff pro- 

cedure, corresponding to c = e~-^/^, yieldsS, Ci = ^^j^. 

We have verifiedi^ that the same result holds within al- 
ternate regularization procedures that do not use sharp 
cutoffs, as long as the Ward identity is satisfied. For 
example, replacing the Coulomb interaction V{r) — > 



with To a length scale and 77 > (putting 



the physical system, at d = 2, sHghtly below its own up- 
per critical dimension) regulates all integrals in a way 
consistent with Eq. (fT7| . We obtain (for rj ^ 0) 



(i)iJ| 



Nauo 



BC 



1 , 4u 
lb \ ?/ lotq 



Im/, 



Nauj 



1 , 4v 
^ , - + In 

16 V 77 ClTq 



IE ) , (25) 
19-67r' 
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with the Euler constant. Once again, while the sepa- 
rate contributions diverge with 77 ^ 0, their sum is con- 
vergent and yields the coefRcient Ci — -'^j^- 

As we have discussed, Eq. ^ impHes that the correc- 
tion to cr'"' is small at low photon energies uj ^ D regard- 
less of the value of Ci. However, at larger (i.e., opticaP'^) 
frequencies, the second term may become significant, de- 
pending on the value of the number Ci . Nair et al? find 
the conductivity to be ct/ct^") = (1.01 ± 0.04). If we 
take the value Ci ~ 0.513 from Ref. ii, however, Eq. ^ 
predicts a large frequency-dependent correction to the 
conductivity that is not consistent with these error bars, 
giving, for example at photon wavelength A = BOOnm (or 
huj = 2.07eV), (7/cr(°) ~ 1.667, assuming the bandwidth 



D = 7.24eV. In contrast, using Ci = 0.0125 of Ref. [lOl, 
yields for the same parameters, cr/cr(°) ~ 1.016, consis- 
tent with the error bars of Nair et al~. In Fig. [T] we 
show the optical transparency that result from both val- 
ues for Ci (using bandwidth D — 7.24eV), along with the 
free-Dirac fermion result as function of wavelength A in 
comparison with experiment. 

Given the smallness of interaction corrections, with 
fj/0-(o) ^ 1.016 for the correct value of Ci (a correction 
comparable to corrections due to the true tight-binding 
band structure) , it is unlikely that optical measurements 
will reveal interaction effects. Electron-electron interac- 
tions are much more visible in the enhanced diamagnetic 
response^ or the hydrodynamic transporlji^. 

In summary, we determined the leading corrections to 
the optical conductivity and transparency of graphene 
and find that they are very small and determined by 
the fine structure constant q;qed up to corrections of or- 
der 1-2% in the conductivity and 0.03 - 0.04% in the 
transparency. Correctly regularizing the UV-divergent 
contributions required using Ward identity arguments to 
resolve previous discrepancies in recent literature (a con- 
troversy that persists^i). Our work demonstrates that 
there are no discrepancies in cr(w) obtained by different 
theoretical methodsU-, such as the Kubo formula, the 
density polarization approach or kinetic approaches, if 
charge conservation is guaranteed at all stages of the cal- 
culation. Our methods confirm and, more importantly, 
justify the result first obtained by Mishchenkoi^ and pro- 
vide a general prescription for calculating interaction cor- 
rections in graphene. 
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